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ANALYTICAL METHOD OF APPROXIMATING THE MOTION OF A SPINNING 
VEHICLE WITH VARIABLE MASS AND INERTIA PROPERTIES ACTED 

UPON BY SEVERAL DISTURBING PARAMETERS 

Bv .James J. Buglia, George K. Young, Jesse D. 

Timmons, and Helen S. Brink worth 


SUMMARY 

An analytical method is presented which approxi¬ 
mates the flight-path deviation of a variable-mass, 
ran able-inertia, spin-stabilized vehicle under the 
influence of ( a) initial pitching motions, (b) thrust 
ro isalinernents, and (c) principal-axis misalinemeni , 
or dy nam ic unbalan ce . I h is method includes a 
first-order analysis of the effect of jet damping on 
the resulting flight path. 

An extensive analysis was carried out on an IBM 
70j electronic data processing machine , by integrat¬ 
ing numerically the six-degree-of-freedom equations 
of motion of a thrusting body with variable mass and 
inertia. This computer analysis was intended to act 
as a check on the method developed herein, and also 
to present other investigators with exact numerical 
answers with which to compare their own results. 
The results from the analytical solution are com¬ 
pared with the results from the numerical integration, 
and the results are found to agree well in all cases 
where no large change in inertia ratio was involved. 

The results are presented in terms of the angular 
deviation in sjuice of the velocity vector from its 
undisturbed orientation , as a function of spin rate, 
for each of the disturbing parameters. It is believed 
that this velocity-vector deviation as a function of 
sjnn rate is a much more intuitive indication of the 
dispersion than the body attitude in space, which is 
the most frequently used parameter in spin-stabdiza- 
tion discussions. 

INTRODUCTION 

With the fairly recent ail vent of very-high- 
altitude rocket vehicles, the engineer has become 
facial with the difficult problem of adequately 


stabilizing his vehicle in the region where the 
atmosphere is so thin that conventional aerody¬ 
namic controls are rendered useless. One possible 
means of stabilization is the use of reaction con¬ 
trols in place of the aerodynamic controls. How¬ 
ever, in view of the diminutive size of some of the 
present-day probes and satellite injection stages, 
a control system of this type would impose a 
severe weight penalty on this final stage, lienee, 
some other means of pitch and yaw attitude 
stabilization must be considered. 

Ballisticians have long known that a spinning 

shell deviates from a predetermined flight path 

much less than a nonspinning shell. It thus 

seems natural to attempt to attitude stabilize a 

rocket vehicle by spinning it and, indeed, this 

method has boon successfully used already in the 

« « 

Vanguard, Kxplorer, and Pioneer vehicles, and on 
Scout and other vehicles currently being designed. 

K' I? n 

Spin, however, is not a cure-all: spinning a vehicle 
gives it a constant attitude in inertial space 4 , but 
not stability along the constantly changing tan¬ 
gent to tiie flight path. Hence, if a large change 
in flight-path angle is involved, the method of 
spin stabilization may not be applicable. Spin 
stabilization is very effective, however, in the case 
of a satellite launcher where the last stage, the 
injection stage, is final at nearly horizontal flight- 
path and attitude angles and it is desired to burn 
out this stage in such a manner that the final 
velocity vector will be nearly horizontal. 

Some investigators in this field (for example, 

refs. 1 and 2) have chosen as a criterion for spin 

stabilization the attitude of the body, defined bv 

* • •• 

tli(‘ two Puler angles in pitch and yaw. It- is fell, 
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however, that a knowledge of the vehicle attitude 
alone is not the best criterion for determining the 
effects of various extraneous forces, because the 
attitude alone does not indicate the final orien¬ 
tation of the velocity vector at stage burnout. 
This velocity-vector orientation is believed to be 
a better indication of the dispersion effects sought. 
Consequently, the data presented herein are given 
as deviation of the velocity vector from its undis¬ 
turbed orientation as a function of spin rate for 
several external disturbing forces. 

The two methods in most general use for deter¬ 
mining dispersion ot a spinning vehicle, with the 
body attitude as a criterion, are the met finds of 
Nicola ides (ref. 1) and Jarmolow (ref. 2). Xico- 
laides determines the dispersion in terms ot angle 4 
of attack a and angle of sideslip & bis parameter 
is essentially ^a*+0*. He also includes aerody¬ 
namic effects, which are neglected here. However, 
he considers the vehicle mass and moments of 
inertia, as well as the applied moment and forward 
velocity, to be constants. His results, therefore, 
cannot be used to study the dispersion of a spin¬ 
ning rocket during its thrusting phase. 

Jarmolow, on the other hand, considers the case 
of variable mass, moments ot inertia, and applied 
moment ; but he 1 discusses only the body attitude 1 . 
Also, analytical results are given only for the pitch 
and vaw rates in a body-axis system. Numeiical 
integration must stdl be list'd to find tbe pitch and 
yaw Euler angles. 

The method presented herein utilizes exponen¬ 
tial approximations for the variation ot the ratio 
of the applied moment- to ihe pitch moment of 
inertia, and for the ratio of thrust to linear mo¬ 
mentum, which appears as a parameter in the 
translation equations. 1 hese approximations 
make possible analytical solutions for pitch and 
yaw Euler angles, and also expressions defining the 
angular deviations of the velocity \eetor. I hese 
expressions include linear damping in pitch and 
vaw. The following additional assumptions are 
made: the ratio of roll inertia to pitch inertia is 
constant: the vehicle is rigid and has rotational 
svmmetrv; then' is no roll moment; only small 
angles are considered; the motion takes place in a 
vacuum. This study was initiated to determine 
the effect of various disturbing parameters on the 
flight path of a solid-propellant, rocket motor 
similar to the final stage of the Scout. 
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SYMBOLS 

constants 

pitch and yaw damping coefficient, 
ft-lb-sec 


constants 

angular-momentum vector of ve¬ 
hicle with respect to center of 
gravity in x b ,y b ,z b coordinate 
system, slug-ftySec 
angular-momentum vector of ve¬ 
hicle with respect to center of 
gravity in x n y c ,z c coordinate sys¬ 
tem, slug-ft “/sec 

angular-momentum vector of ve- 
hide and jet with respect to 
center of gravity in x n y n z r coor¬ 
dinate system, slug-ft'/see 
unit vectors along the A*-, 3 
and Zf-axis, respectively 
unit vectors along the AV, 3,-, 
and Z f -axis, respectively 
V('hi<*li‘ moiiH'iits of int'rtia about 

AV, 3 V, and Z,-axis, respee- 
tively, slug-ft- 

J y or 1 A (where slug-ft J 

vehicle products of inertia about 
AV, 3V, and //-axis, respec¬ 
tively, slug-ft" 
constants 

length from vehicle center of grav- 
ity to nozzle exit plane', ft 
disturbing-moment vector, J/ A i c 

-r AfyjA H-ll> 

magnitude of disturbing monu'nt, 

ft-lb 

mass of vehicle, slugs 
element of mass, slugs 


M„ Jl at. /—<>, I/sac 2 

/.v 

rad tans/sec 


/V/Z 


angular-velocity components of 
vehicle about AV, IV, and Zf- 
axis, respectively, radians/sec 
unless otherwise indicated 
distance from center of gravity to 
jet exit, /v.a'L-H - , yjrT>Y.;;k r , it 
magnif ude of r,, (1 
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t time, see 

t hrust, lb 

vehicle velocity vector, ft /sec 
magnit tide of V, ft/sec 
jet (exhaust velocity, h". <A -i r + 

Ve.rJe+Vr.zK, f t/SCC 

magnitude of V 0 ft/see 
coordinate's measured along hodv 
AV, ) V, and Z,-axis, respec¬ 
tively 

X f)} ) tiy Zb bodv-tixe'd system of orthogonal 

axes, origin at body center of 
gravity 

j' n ij n z f coordinate's measured along AV, 

>V, and Z ( .-axis, respectively 

X r ,Y n Z c systemi of orthogonal axes which 

pitches and yaws with vehicle 
but does not roll, origin at body 
(•enter of gravity 

coordinates of (‘lenient of mass 
causing dvnarnic unbalance 

p 4 

angle' of attack, radians 
angle of sideslip, radians 
orientation of velocity vector in 
vertical plane', radians 

6 thrust inisalinement angle', radians 

(except as noteel) 

e principal-axis misalinomont angle, 

radians (exempt as noted) 

X orientation of velocity vector in 

horizontal plain 1 , radians 

a total oriental ion of velocity vector, 

radians (except as noted) 

0 pitch Euler angle*, radians (except 

as noted} 

yaw Euler angle*, radians (except 
as noted) 

angular velocity of nonrolling svs- 
tern in inertial space, 
o) Xc kn radians/si'c 

o)y ( .,c components of <a r along )V and 

Z,.-axis, re*spe*ct ively, radians/sec 

Subscripts: 

o denote's ([uantity at zero time 

A,) ,Z ronipoiinits parallel to AV, 3 V, 

and Z A -axis, respectively 

AV, )V,Zc eompoiH'iits parallel to AV, I V, 

and Z f -axis, re'spectively 

A dot over a quantity demotes (Inferentialion 

with respect to time. 




CC 


ANALYSIS 

The mat hemal ieal analysis of a spinning rocket 
differs in several respects from that of a spinning 
projectile. Jn the first place*, the rocked lias mass 
and inertia properties which vary with time, 
whereas those* of a spinning proje*e*tile are constant. 
Secondlv, the rocket can have a moment dm* to a 
thrust inisalinement. This moment, which is 
fixed to the spinning vehicle, give's rise* to motions 
not found in the* analysis of projectile motion. 
The motions due to an initial pitch or yaw rate 
and a dynamic unbalance* an* common to both 
type's of boelie'S, but the analysis becomes more 
involve'd when a variable-property vehicle is elealt 
wit h. 

The e*quations of motion programed on the 
IBM 704 e'lectronie data processing machine* for 
comparison purpose's in t his si inly are* t he* standard 
six-di'gree-of-freedom equations found in most 
aeronautical textbooks (for example, ref. R). 
The ('([nations programed, along with the* various 
transformations that were used, are presented for 
reference in the* appendix. 

An approximate analytical method of predicting 

the deviation of the velocity vector dm* to the* 

* 

effects of thrust inisalinement, initial pitch rate, 
and principal-axis misalinement will be* developed 
in this section. Two coordinate* systems will be 
used in this analysis. The* first ,:y) is a 

body-fixed coordinate system with axes fixed at 
the center of gravity of tin* vehicle* and the* 
coordinate* x f) extending along t la* body longit udinal 
axis. The second coordinate system (./■,. ,y n .z c ), in 
which the* equations of motion will be* derived, is 
referred to a body-fixed system of axes which can 
pitch and yaw wit 1 i the* vehicle* but does not roll 
wit hit. The t wo sots of coordinate* axes arc* a lim'd 
at /.=_(). (See* fig. 1 .) Pile* pitch Euler angle* 0 
and the* yaw Eule*r angle ^ are assume'd to be* small 
so t hat the order e>f the Euler rotation is immate*- 
rial. This assumption of small angles simplifie*s 
the* integration for 7 and A, the flight-path angle's 
in the* vertical and horizontal plane's, re*spe*ctively, 
and also permits these* flight-path angles to 
combine* vectorially to give* a, the* angle between 
the* velocity vendor and some* fixed reference line, 
take'll here* as the* longitudinal axis of the vehicle* 
at zero time. 

The* angular-momentum vector of the* vehicle 
in the x b ,}f b) z b system instantaneously fixe*d in 
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Fidi'KK 1-Coordinate system and definition of symbols used in analytical development. 


inertial spare can be written 

H t iJ A p f\ y(J ~T (4j q I yyj' I\ y}>) j ft 

IIYz f j)k i> ( 1 ) 


Xo\\\ in the non rolling system (x n y c ,z c ), which 
has anyular velocity in inertial space, 

CO ( — (0, CO Ycy M'S. c) — W V (- j c ~l' UzMc (" ) 

The conation of motion about the vehicle center 
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of gravity is 


(I Hr 
(It 


M 


(:*) 


The rate of change of angular momentum with 
respect to the vehicle center of gravity including 
jet effect is 

(/Hr (IH r . . 

~l j~ = ^ **)“ « c X H r + m r e X (w <- X r f — \ ,.) (4) 

where w f XH ( . takes into account tlie* rotation 
of the coordinate* system. 

it is usually conventional in missile work to 
choose the body axes as principal axes and to 
assume a symmetrical vehicle, so that 


Iy-Iyr I 
Ix y — I yz — J xx ~~ 11 


(o') 


With the •se assumptions, equation (1) may lx* 
written in the nonrolling system as 


H f - TxpK + I<*> yc } e -\ t^zMr 


( 0 ) 


If the* further assumptions are made that the 
jet is located on the AVaxis and the jet exit 
velocity relative to tlie* nozzle exit is directed 
along t he* negative* AVaxis, then 


and 


1 f, AV "" 1 t 

1 (. Yc~~ 1 r.Xc~~ d 
G. xc~ if 

r f Yc~ r <\Zc~~ fl 


(7) 


( 8 ) 


Expansion of equation (d), with the aid of sub¬ 
stitutions from equations (2), (4), and (0), yields 

( I:\P ~U / xp) b "b ( /u> Y c -r 7a> y P +1 x p P -p th rj cc Vc ) j r 

r IIxpv>Yc~\- )k c —- IV1 (9) 

If damping in roll is neglected (that is, if / x -- 0), 
equation (9) becomes, in component form, 


IxP-Mx 

1°- Yc^~ JxP^zr^C LOy c - - Af y 

l^Zc 7 v P^y f -t" C iCy c — A f y 


(10a) 

(10b) 

(10c) 


whe 


re 


c'= mr ;;+/ 


The quantity c'fl( — c) will be considered constant 
for this analvsis. 

i 

If tlie* assumptions are made that the disturbing 
moment lies in the vertical plane* at ze*ro time 
and that the*re* is no moment about the* roll axis, 
t hem 


yfx 

My 

M x 


= 0 


Af cos pt 
Af sin pt J 


( 11 ) 


and equations (10) become 

Jxfi — 0 


TuyrA" Ix/>Uzc + (-'uy c — M COS pt 


(12a) 

(12b) 


Iuzc — IxpwycAc'uzc~ Af sin jd (12r) 
Eejuation (12a) immediately shows that 


p- ( V)nst ant jt 0 


Let 


h 

I 


Pv=H 


(l :>) 


and 


\f ( M 

j 1 i> f 


(14) 


It is assumed that the ratio [ x jl is constant; Af 0 
and k are* also constants, found by fitting an 
exponential to the function Afjf as shown in 
equation (14). r Plu*ii equations (12b) and (12c) 
become 

a >y e -\- no) X f-rCa>r c ^ Af a e kt cos pj (Ida) 
Vc — nuve -r cu Zc — AI y i sin pj (ldb) 


'These* equations can be integrated simply, and 
with the initial conditions «»(())-■ a' rr .o, a? Zr (d)=0 
give 

(cijjv „— ( 1 i )e~ ct cos nt — (\e~ ct sin ni 

cos pJ-r(\e Kt sin pj (lCni) 


cczc~(^vc,u — W lt sin nt 1 ( X ct cos id 

- XX' sin p„t—(\e kt cos pj 


where* 




Af ti (Po — n) 

(c + k)-± (po—n) 2 


(10b) 


(17a) 
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U,(c + k) 


r> 


(c+t) 2 +(y ^“ w ) 2 


(17b) 


Now, bv virtue of the original assumption that 
the pitch and yaw motions are small, and as a 
consequence of the coordinate system used, the 
following expressions can be written: 


0— fu>Yc<M 
f UZftH 


(ISa) 

( 1 *b) 


Substitution of equations (10a) and ( 10 b) into 
these integrals with the initial conditions 
0 ( 0 ) ™^( 0 ) — () gives 

e—A(e kt cos pj— \)-\-Iie kt sin pj 

— l)(e~ ct cos /// — 1 ) -\-Ke~ ct sin id (19a) 

\p~ Ac*' sin pj — B(e k( cos pj— T) 

— })(~ c: sin id — FA(~ ct (‘os nt — 1 ) (19b) 


w 


here 


.1 


II 


( ,k'-_ ( '\P- 

F \ pi 

< '« Pop'll' 

F-{-pl 


{20a 


n 


eUOjv.'/' 2 ) a ( < 


r- 4- “ 


(20c) 


PF- 


idjjL'y,. ,, - { y 1 c( 1 

c 2 -[-/c 


These equations give time histories of tlie body 

attitude angles in an inertial coordinate system 

taken to be coincident with the body coordinate 

system at zero time. These angles 0 and \p are 

the' angles commonlv used to discuss the effects of 

spin on the stability or dispersion of a vehicle. 

However, body attitude' alone does not adequately 

describe the complete motion in space. The act ual 

dispersion of the vehicle, defined herein as the 

angle between the velocity vector at zero time and 

the velocity vector at burnout, is a much more 
% 

meaningful parameter because it depends on the 
translation of the vehicle', and lienee on the quan¬ 
tities which define its translation history, as well 
as on the rotational parameters. Tin* preceding 
argument can be substantiated from physical 


considerations. The equations defining the rota¬ 
tional motion of the vehicle are completely inde¬ 
pendent of the translational motion. (See, for 
example, eqs. (Al) to (AG).) Hence, for given 
inertia characteristics, spin rate, and disturbance, 
the rotational motion of the vehicle is completely 
defined. Now, in order to get the actual transla¬ 
tion in space, the body attitude parameters must 
be combined with the translation parameters; 
for example, equations relating quantities such as 
thrust and mass ratio are integrated to give the 
translational motion. As an illustration, it is 
certainly possible to have a spinning vehicle point¬ 
ing vertically upward and moving horizontally at 
ignition. Then as the motor fires the velocity 
vector climbs above the horizontal. Thus, the 
final direction of the velocity vector has changed— 
the dispersion as defined here is reduced—hut the 
body attitude has remained constant. An effect 
such as this one is not evident when only the rota¬ 
tional motion is defined. 

Jn order to determine 7 , the orientation of tin* 
velocity vector in the vertical plane, reference is 
made to tin* normal equation of motion of figure' 
2(a). Since a — 6 — 7 , if only small angles are 
taken ii to account this equation becomes 


<h , r 

,it + n,V y 


=dde- 


T 


0 cos /V, 


( 21 ) 


mV v ; m\ 

In a sii filar manner a differential equation for X, 



(a) 

mVy - Sin 8 cos p Q t COS a + T cos 8 sin a 
mVy = TS cos p Q t + 7"(0-y) 

(a) Vertical plain*. 

IToi rk 2. Equal ions of motion 
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the' orientation of the velocity vector in the hori- 

« 

zontal plane, can be obtained from figure 2(b): 


d\ T 
dr mV 


A 


T T 

mV^^mV SS1 " vJ 



Since only the changes i11 7 due to external dis¬ 
turbances rather than (lie natural change in flight 
path due to gravity are of interest, the effect of 
gravit v is omit t ed. 

Now, in order to integrate equations (21) and 
(22), the quantity TjmV must be put into a form 
that is easilv integrated. It can be shown that a 
plot of TjmV iVA a function of time can be approxi¬ 
mated to a sufficient, degree of accuracy bv an 
exponential function. 'I'lms, define. 


and 


T 

mV 








mV\- T sin 8 sin p Q t cos & + T cos 8 sin 
m V X - T 8 sin p Q t + T )8 


In the process of integrating equations (21) and 

(22), the quantity VA(M k2t appears in several in¬ 
stances. Therefore, to further simplify Ihe inte¬ 
gration , it was assumed that 

~c ( 24 ) 


(h) Horizontal plane. 

14 <• u kh 2 . Concluded. 

If equation ( 2 . 4 ) is substituted into equations 
(21) and (22) along with (lie expressions for 6 and 
\p from equations ( 19 ), integration of the resulting 
linear differential equations, utilizing equat ion ( 24 ) 


and tlie initial conditions 7(0) A(0) 0, yields 


y — ( 

r/ ' 

~Ae «•*+*''((',+/-) - p„Be w - «»* Se c '*‘( 

(r y k y+p i 'Pi+p l- 

eos p,J + 

Ap„e u 


k) 


5 e Vs ‘p- 

n+r 

. r/ >(('■, — c)e ,( 'i~ cU -\-Fne 
sin p„t—\ ' - — --• - •• 

L (fv-c)-+«- J 

cos ni — 

'nDe"'-~ 

c "-Fe , ~‘- r,, (/--C) 

(F,-c)-+n? 

<>» 

i—i 



\-k)-Bp a 

TV J L (p 2 +W+Pl 


I)((' 2 -c)+En_A 1) d( f 2 

(Tt-c) 2 +n* V 2 r 2 a +;rj / 



- k. 

\=<>r ke 


U __ ( 

V('A 


— p„Ae «*+*>«_ B{(\-\-k)e t( ~*+ ku pM v *‘ 


(<V\ ky+ri 


‘ ,,p 


U^+ k >'((’, + k) - p„Jie ,Vs+ * H 

{V+kyipi 



cos nl 



eU (r,- o+uFa'" ! 
(P-c)--i 


sin nt 



Ffvi _r- p,A-J{« 1 k) 

r, J L i /•)= | pi 



nl)-(r,-c)F n FA I 

e) 2 - 1 ,- it- + 7 y 7v_l/ 



021963 —02 
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Thru, because of llie assumption of small angles, 
7 and X ran la* combined vectorially to give <r , ihe 
angle between tin* fi.\(‘(l reference line and the 
velocity vector: 

<r-y?+X 2 (27) 


with the assumption c<< 1, substitution of equa¬ 
tion (21) into (‘([nation (20) gives 



A disturbing parameter consisting of a princi¬ 
pal-axis misalinement can la* bandied in the fol¬ 
lowing manner. Assume that the mass element 
causing the dynamic unbalance can be placed oil 
the center line of the vehicle a distance t/ u from 
the center of rotation, and that it lias a mass 
m u , as in sketch l. The moment about the renter 


Yb 



A/. vU tt y u m n (2S) 

Now, the produc t of inertia due to this mass is 

Lxy XuVuMu ( 2 b) 

Hence, the moment can be written in terms of the 
product of inertia: 

\f pul.xy (dO) 


It may be desirable to (express the product of 
inertia in terms of the principal-axis misaline- 
ment angle c. From the properties of moments 
of inertia it is known that the angle e between 
the principal axis and the axis about which the 
moments and products of inertia are computed 


is given bv 


tan 2< 


2/.vr 

Ix-I 


( 21 ) 


If, as will be the usual case, the product-of-inertia 
term is due to construction technique or placing 
of instrumentation, the product of inertia will be 
constant, the small change due to the center-of- 
gravity shift during burning being neglected. In 
this case, e will be a function of time. The time 
history of e can be computed from equation (21) 
since / v and / are known at anv time 4 . Then 

' v i 


If this expression is used in (‘([nation (14), and M (t 
and k are found by a curve lit, solution of capta¬ 
tions (25) and ( 2 b) will yield histories of 7 and X 
due to a principal-axis misalinement e. 

It should be noted that, to the degree of approxi¬ 
mation used herein, the effects of the various 
disturbing parameters can be combined linearly 
to yield the total dispersion. For example, equa¬ 
tions (25) and (2b) may be written as a linear 
combination of terms each of which is a function 
of ojyc ,1 , A/«, or 8 but not of any t wo of the variables. 

It should be further noted that, again to the 
degree of approximation assumed here, the angles 
7 and X, and consequently the total dispersion <r, 
vary linearly with a ) YCt0 and Aif taken sepa¬ 
rately, so that it is not necessary to make a 
separate calculation for each value* of moment or 
pitch rate desired. This linearity is supported 
by the calculations performed on llu* IBM 704 
compuler, as illustrated in several of flu* figure's. 

RESULTS AND DISCUSSION 


A comparison of the results of the numerical 
integration (on the IBM 704 computer) of the 
exact, -‘(plations of motion given in the appendix 
with tlie results of the analvtical solutions do- 

t ■ 

ve loped herein is made in figures 2 to Hi. The 
quant-ivies 6 , \p, 7 , X, and a were calculated from 
(‘([nations (M)a), (10b), (25), (2b), and (27), 

respectively. The constants used arc given in 
tables I and II. Figures 7 to lb are plotted for 
f = 20 seconds (burnout). These figure's serve a 
double purpose. Since not everyone lias access 
to an IBM 704 or to a computer of sufficient 
capability and storage capacity to permit the 
progra ning of a problem as complex as a six- 
degree-of-freedom rigid-body analysis, many in¬ 
vest ig» tors have probably developed analytical 
results without being able to check the accuracy of 
their methods. Thus, the figures presented can 
serve to check these methods and indicate regions 


of validity and applicability. At the same time, 
the figure's serve* to illustrate the accuracy of the 
analytical method developed herein. As can In* 
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TABLE I. PH YSICAL AND PER Ft >R M ANCE 
(41A RA(TE R1STI(• K OK ROCKET MOTOR 
rSKI) IN CALCULATIONS 


Hull moment of inertia, slug-ft 2 
Pitch (yaw) moment of inertia, slug-ft 2 

Mass, slugs ....... . 

Applied moment, ft-lh (5 in radians) 
Thrust, lb ... . 

Burning time, see 

Init ial velocity, ft/sec . 

Final velocity, ft/sec _ 


. 5.25 0.18/ 
- 84.15 0.88/ 

18.05 0.18/ 
10,7105 • 1505/ 
8,000 
50 

10,000 

10,810 


('output cd constants (for 0.10° t lirust misalinement): 


_ . . 

X. .. . 

(\ _ 

r 2 _ 

5, radians ..... _ 

to» radians/sec 
Ix'I 

1/sec 2 .. 

/;. ..._ ... 


0.01585 

0.0200 

0.008 

1.880 

0.050 

0.001745 

0 

0.104 

0.5470 

0.0502 


TABLE II. PHYSICAL AND PERFORMANCE 

characteristics of spherical rocket 

MOTOR 


Roll moment of inertia, slug-ft 2 

888.ti 10.0/ 

Pilch (yaw) moment of inertia, slug-ft' 2 

408.0 10.0/ 

Mass, slugs . 

10.88 0.51/ 

Thrust, 11) 

8.400 

Burning t ime, see 

80 

Moment, ft-lb 

5 | 0.01/ 

Initial velocity, ft /see 

10,000 

Final velocity, ft,/sec 

20,050 

Computed constants: 


^*i -----. 

0.02125 

k 2 _ -. - 

... 0.0700 

A*. .. ... . 

0.200 

£i _ 

1.280 

C 2 ___ . .. 

.... 0.107 

5, radians ........ .. . . 

0 

(i)Yr o, radians,/sec .... 

0 

Ix jl ____ . 

0.08 

M n , 1/sec 2 __ _. . 

_ 0.0124 

k --- -- 

0.042 

seen from a general perusal of figures I 

\ to 15, the 

accuracy obtained is quite good in 

the range 


invest igated. 

Figure 3 shows the pitch Euler angle 0 plotted 
against the yaw Euler angle i£, for the last few 
seconds of burning of the rocket motor, under tin 1 
action of a moment, produced by a thrust misaline- 
ment 5 of 0.10°. This figure shows, essentially, 
the attitude of the body in space during this time. 
If ran be seen that, although the analytical method 
predicts frequencies of nutation and precession 
which are slightly high, the amplitudes of both 



^,deg 


Florin; 8. Undamped body attitude. />„—80 rpm; 

<Y~ 0.10°. 

motions, which are 1 Ik* parameters of primary 

interest, are- predicted rather well. This curve 

does not include the effect of jet damping. 

Figure 4 shows lime histories of tin 4 deviation 

of the velocity vector a. The analytical method 

« * 

predicts results which are slightly conservative. 

In figures 5 and (\ arc plotted time histories of 
the pitch Euler angle 6 and the yaw Euler angle \j/ y 
respectively. Here if will he noted that, while 



Fiiaiu; 4. -Deviation of velocity vector as a function of 
time, for a thrust misalinement angle 6 of 0.10°. Un¬ 
damped body; p 0 — 80 rpm. 
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Deviation of velocity vector, <r, de 






Deviation of velocity vector, cr, deg 
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Figure 


8.—Effect of spin rate on the deviation of the velocity vector due to a thrust 

Damped body; t ~30 seconds (b irnout). 


niisalineinent angle* 5 of 0.10°. 



Deviation of velocity vector, <r,deg 
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Khjuuk U Linear t;fTcu*t of thrust misaliuement on the 
deviation of tin* velocity vector for several spin rales. 
Undamped body; /- 30 seconds (burnout). 



Fine itK 10. 


Fifed 


of spin rate 
angle 


on the deviation of the velocity vector due to a principal-axis 
e of 0.10°. Undamped body; ( lit) seconds (burnout). 


misalinenn 





Deviation of velocity vector, a, deg 
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1t>. -Effect of spin 


rate on the deviation of the velocity vector due to an applied moment for a 25-inch spher¬ 
ical rocket motor. Undamped body; C--50 seconds (burnout). 


CONCLUDING REMARKS 


The of roots of spin stabilization on the flight - 
path angular deviation caused by initial pitching 
ratty thrust- misalinement, and principal-axis inis- 
alinenient have been evaluated by an approximate 
analytical solution which yields the flight-path 
angular deviation as a function of spin rate for a 
vehicle whoso mass and inertia properties are 
varying with time. Agreement with a numerical 
integration of the six-degree-of-freedom equations 


of motion was found to be good in all cases in 

which the inertia ratio was essentiallv constant. 

* 

Spin stabilization was found to he very effective 
in averaging out the errors caused by initial pitch 
rules and thrust misalinemonts, but was com¬ 
pletely ineffective in controlling the errors in flight- 
path angle due to principal-axis misalincnient. 

Langley Research Center, 

National Aehonattics and Space Administration, 
Langley Am Forte Base, Va., April Jl, UWl . 






APPENDIX 

EQUATIONS OF MOTION IN BODY-FIXED AXIS SYSTEM 


The symbols used in this appendix are defined 
as follows: 




Ay. IY'lz 


Ay z 

I 

i iy D?' tl i 

\I x MyM z 


m 

]) y qj' 

r et X 

[ A, 


V 

xx y. 


xr,z 


7 

0 


X 


a 


components of external force along 
the body A'-, E-, and Z-axis, 
respectively, 11) 

moments of inertia about body 
A'-, Y- y and Z-axis, respectively, 
slug-ft 2 

product of inertia, slug-ft 2 
I y or ] y (where I r —I x ), slug-ft 2 
direction cosines between inertial 

axes and body-fixed axes, v— 1, 

• * ‘ 


external jnoments about body A-, 
Y- % and Z-axis, respectively, 
ft-lb 

mass of vehicle, slugs 
angular-velocity components ol 
vehicle about body A-, )and 
Z-axis, respectively, radians/sec 
longitudinal distance 4 frojn vehicle 
center of gravity to jet exit, ft 


t runs form at ion mat rices 
velocity components of vehicle 
along A'-, and Z-axis, respec¬ 
tively, ft see 

vehicle velocitv vector, ft /sec 


orthogonal bodv-axis svstem 
coordinates of (tenter of gravity of 
vehicle along A-, ) and Z-axis, 
respec-tively, ft 

fixed orthogonal axis system, coin¬ 
cident with bodv-axis svstem at. 

• 4 

zero 1 ime 

orientation of velocitv vector in 

* 


vertical plane, radians 
pitch Killer angle, radians (except 
as noted) 

orientation of velocity vector in 
horizontal plane, radians 
total orientation of velocity vector, 
radians (except as noted) 


<t> roll Euler angle, radians (except as 

noted) 

\j/ yaw Euler angle, radians (except 

as noted) 

A dot over a quantity denotes differentiation 
with respect to time. 

The rigid-body equations of motion as pro¬ 
gramed on the IBM 704 electronic computer are 
the standard six-degroc-of-freedom equations of 
motion found in most aeronautical textbooks on 
dynamic, stability (for example, ref. d, eqs. (10-11)), 
with damping in pitch and yaw added to the 
second and third moment equations, respectively. 
These 1 equations are repeated here for reference. 
The tluec force (‘(plations arc 4 

m (it ~rr-\ w q) A ( A1) 

m (v—wp + nr) — X F y (A2) 

m (w— u <i + rp) — XI Vz (Ad) 

and tin moment equations including jet damping 
are 

/ •■/) Ix/j \~ Hz Iy) ( P‘ Iy? !t<i~ d Ay (A4; 

/ if/ j (/ I z )rp-y Ay z(J> 2 — >' 2 ) 

m 

4-/y(/-f fhr; v(/- Xd/r (Af>) 

fyj' Ay 'P + (I r Ay ) P ( l A Ixz ( p‘ 

* Iyf-\ thrift--- 'ZZMy (AO) 


The plant it ies us<‘d liere ai*c 4 r( 4 f( 4 n*( 4 d to a set of 
axes ib ed to and rolling with the bodv, whereas 
the analysis in the body of this paper is referred 
to a set of non rolling axes. 

The quantities u , v\ w and />, q % r are found by 
iutegraion of these 4 equations by tin 4 Kunge- 
Kutta nethod, and are then transformed to a set 
of iner ial axes. The following assumptions are 
made t > simplify tin 4 integration: 

M >ments and products of inertia vary 
inearlv with tinn 4 

4 


IS 
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Mass varies linearlv with time 

* 

o V 

Damping coefficients (m/j x ~\-1) are constant 

The inertial quantities are found by means of 
the following transformation equations (ref. 4): 


The angular rates */, r about the body axes 
are transformed to angular rates 0 , 0 , \p about the 
space axes by means of the following transforma¬ 
tions (ref. 4): 


.r 


7 / 


0 =p-|-tan 0(7/ sin 07 r cos 0 ) 

// 

=r/ , iMr,,][7’ :i i 

r 

(A7) 

0 — q cos 0 —/* sin 0 ^ 

1 

-- (7 sin <p-\-r cos 0 ) 

COS 0 

— ■ — 





(A 13) 


whore the transformation matrices [7 T ,], [ 7b], and 
[7a] are given by 


[7i] 


cos 1 p —sin 1 p 0 

sin \p cos \p 0 

0 0 I 


(AS) 


i n 


cos 0 0 

0 1 
-sin 0 0 


i r,\ 


1 


sin 0 

0 

cos 0 _ 
0 

sin 0 


(Ah) 


0 cos 0 

j_() sin 0 cos 0 J 
Multijjl\'ing the matrices gives equation (AT) as 


(Alt)) 


x 

* 

Ti 


/1 M I ii 1 


1 2 iff 2 M-2 

Ua "b H*. 


U 

r 

V' 


A1 1) 


where 


/. — cos yp cos 0 
/ 2 —sin \p cos 0 

! :i — — sin 0 

/ft T ~ cos yp sin 0 sin 0 —sin \p cos 0 
/k 2 — sin \p sin 0 sin <j>-\ cos \p cos 0 Y (A 12 ) 
//<3 = cos 0 sin 0 

//j = cos \p sin 0 cos 0 -{ sin \p sin 0 




Appli cation of these transformations to the body- 
axis quantities ?/, /\ >c and p, (/, r yields time 
histories of 7, 7/, ” and 0 , 0 , 1 p. Then, from sketch 
2, 7 , X, and a can be found: 





Sketch 2. 


Thus, 


7 — tan 1 


i 


\x 

7/ 

7 


A —tan 


-1 


4 _ 1 ’V' 7^ 


<7 — t a 11 


x 


/E> = sin \p sin 0 cos 0 — cos xp sin 0 
M j —— (‘OS 0 COS 0 
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